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ABSTRACT 

We have used deep /-band (F814W) images from the HST archive to study 
the globular cluster systems around the brightest cluster galaxies (BCGs) in 
Abell 262, 3560, 3565 and 3742. Three of these BCGs have inner dust lanes and 
peculiar structural features which indicate past histories of low-level interaction 
and accretion. The deep /—band WFPC2 images have photometric limits which, 
for all four galaxies, reach near or just beyond the GCLF turnover point. Their 
specific frequencies are 8.24 ± 1.65, 4.66 ± 0.93, 2.58 ± 0.52 and 2.62 ± 0.52 
respectively, all within a factor of two of the normal range for giant ellipticals. 
We obtain new estimates of the GCLF turnover magnitudes, which are shown 
to be consistent with an adopted Hubble constant of Hq ~ 70 km s~^ Mpc~^ on 
the "Hubble diagram" of GCLF turnover apparent magnitude versus redshift, on 
a distance scale where the fundamental GCLF calibrator E galaxies (M87 and 
others) in Virgo are at = 16 Mpc. 

Subject headings: galaxies: elliptical and cD — galaxies: evolution — galaxies: 
individual (A262, A3560, A3565, A3742) — galaxies: star clusters — galaxies: 
structure 



1. Introduction 

Globular clusters are found in their largest numbers within giant elliptical and cD galax- 
ies. While the clusters constitute only a small fraction of a galaxy halo's mass, they are the 



^ Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with proposal ID 5910. 
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most easily identified surviving structures of the initial stages of galaxy formation. This 
makes them important tracers of galaxy formation models. From deep wide-field imaging, 
several features of the globular cluster system can readily be derived, including the radial 
distribution, luminosity distribution, and the total population (specific frequency S'jv). 

In this paper we obtain measurements of the globular cluster systems (GCSs) in the 
central giant E galaxies within four Abell galaxy clusters in the rcdshift range cz ~ 4000 — 
4700 km s-^ These are NGC 708 in A262, NGC 5193 in A3560, IC 4296 in A3565 and NGC 
7014 in A3742. The raw data are WFPC2 /—band exposures drawn from the HST archive. 
Preliminary studies of these systems were made by Laucr ct al. (1998), but using only the 
PCI chip data; we make use of all the WFPC2 data available, more than doubling the 
amount of information on the globular cluster luminosity function (GCLF) in each galaxy, 
and also enabling us to study the GCS radial distributions and specific frequencies with 
much higher confidence. Finally, using the GCLF turnover magnitude and redshift of these 
galaxies and others in the literature, wc construct a purely GCLF-based Hubble diagram 
and show that it is consistent with i/o = 70. 



2. Data Reduction 

The original images were taken in the F814W filter by Lauer et al. (1998) with the 
WFPC2 camera, between January and April 1996 (Program ID 5910). The target galaxies 
were centered in the high resolution PCI chip, hence maximizing the total globular cluster 
population falling within the field of view of the entire camera. The exposures (totaling be- 
tween 9200 and 16500 s) were sub-pixel-shifted. With standard IRAF and STSDAS routines, 
wc CO- added each group to reconstruct clean composite images free of bad pixel artifacts and 
cosmic ray contamination. 

Various small areas on each image were masked out in subsequent analysis to eliminate 
contamination from diffraction spikes of bright stars, or neighboring galaxies in the field. 
Also, the innermost parts of the target galaxies themselves were masked out, within radii of 
3. "7 , 3. "8, 5. "8 and 2. "8 from the centers of the galaxies respectively. 

We next subtracted background light from the combined images. For this purpose we 
used a simple median filter on the WF2,3,4 fields, where the background light gradient is 
low. However, for the PCI chip, we fitted elliptical isophotes (using STSDAS and ELLIPSE) 
to the galaxy, generating synthetic isophotal contours which were then subtracted from the 
original. 

With the galaxy light subtracted, we found that some of the galaxies exhibited various 
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types of interesting substructure. The most prominent of these is the inner region of NGC 
708, which exhibits an easily seen dust lane (see the illustrations of each one shown in Lauer 
et al.). Smaller dust lanes were also detected within NGC 5193 and IC 4296, indicating 
likely accretion events in the recent past. Giovanelli et al. (1982) found that A262, the host 
cluster for NGC 708, contains galaxies deficient of HI gas, suggested as due to ram-pressure 
stripping. This gas could be accreted by the central cD galaxy, with its large dark halo. 
A262 is a spiral-rich cluster, characterized by a X-ray source, 3U 015-1-36, which appears 
concentric with the central galaxy. The gas and galaxy virial temperatures for the cluster 
have recently been shown to be comparable by Neill et al. (2001). There seems to be a 
general pattern for HI deficiency to be associated with X-ray emission, some other examples 
being A2147 in Hercules, or Coma. 

Our procedures for photometry of the four combined images follow the basic sequence 
outlined in more detail elsewhere (e.g., Woodworth & Harris (2000), Kavelaars et al. (2000)). 
Of the many detected faint starlike objects on the frames, the vast majority are globular 
clusters within the BCGs. At these distances (~ 55 Mpc) they appear as unresolved point 
sources, even in the PCI frames. Hence, it is easy to perform conventional point-spread 
function (PSF) photometry on the frames. Independent empirical PSFs were constructed 
from several moderately bright, uncrowded stars in each of the four WFPC2 frames with 
DAOPHOT. Then ALLSTAR (Stetson 1994) was used to generate the final photometry. A 
standard detection threshold of 3.5 times the rms scatter of the sky background was adopted 
for the image detection. Crowding at all magnitudes was completely negligible in these 
high-latitude fields. 

An image classification algorithm (CLASSIFY; defined by Kron (1980) and Harris et 
al. (1991)) was then used to calculate radial image moments of the candidate objects and 
thus to separate stellar from nonstellar objects in an objective way. By using artificial-star 
data passed through exactly the same measurement process, we established boundaries for 
the CLASSIFY parameters such that at least 95% of the true stars were retrieved at all 
magnitudes (for very similar examples with illustrative graphs, see Woodworth & Harris 
(2000) and Kavelaars et al. (2000)). As usual, for the faintest objects the image moments 
become very uncertain, and it becomes difficult to distinguish between stellar and nonstellar 
objects. Also, some nonstellar objects will be accidentally classified as stellar. These are 
statistically removed from the final GCLF by subtracting a background luminosity function 
(see below). 

Finally, to define the photometric completeness function we constructed annular rings 
around the BCGs, and employed extensive artificial-star tests to measure the detection 
completeness /(m,r), as a function of magnitude and radial distance. In practice we found 
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that the PCI chip with its strong background hght gradient was the only area where / was 
a significant function of radius; on the WF chips a single function /(m) could be used. 

The instrumental magnitudes were converted to the Johnson- Cousins / system with the 
standard transformations for F8141^ found in Holtzman et al. (1995). For the four indi- 
vidual BCGs, Galactic extinction corrections of Aj — 0.16, 0.10, 0.11 and 0.06 respectively 
have been adopted {Aj is given by Ai = 1.82 Eb-v, where the Eb-v values are obtained 
for each galaxy from NASA/IPAC Extragalactic Database (NED)). To employ the trans- 
formation equations and also to step back and forth between / and the (more normally 
used) V magnitude scale for globular clusters, we have simply assumed a color index of 
{y — /)o = 1.1 ± 0.1, typical of moderately metal-rich globular clusters in giant E galaxies 
(e.g. Kundu et al. (1999)). The intrinsic range in {V — /)o, folded through the transfor- 
mation equations, will not introduce uncertainties larger than ±0.03 in the calibration of /. 
The assumed {V — I)q value is the mean value representative of most other gE galaxies, and 
if the GCSs were entirely metal-rich or metal-poor, the error introduced by the assumption 
would be at most 0.1 mag. 



3. Analysis and Discussion 

3.1. Radial Profiles 

The projected number density a of detected objects around each galaxy plainly reveals 
an extensive GCS concentrated around the galaxy center in each case. The profile is rea- 
sonably well represented by a simple power- law form a{r) = o"cK^) + '^hg = a + (Tbg, where 
abg is the background number density of starlike objects (mostly faint, small galaxies which 
passed through the image classification routines, plus a few foreground Galactic stars). To 
obtain the profile parameters of the GCS for each galaxy, we subdivided the WFPC2 fields 
into annuli 50 pixels wide, centered on the BCGs. The number density of objects was then 
calculated down to a cutoff magnitude at which incompleteness corrections were still small. 
The projected number density is then just a = N/A, where N is the number of detected 
objects within each annulus and A is the area of that annulus which falls within the WFPC2 
boundaries (minus the small masked-out areas). Completeness corrections, though small, 
were explicitly accounted for. Finally, the background density a^g on each of the four fields 
was defined as the mean of the outermost eight annuli, which fall on the outskirts of the 
WF chips. This corresponds to a radial distance greater than 105" (or about 30 kpc) from 
the centers of the BCGs. Although the GCSs probably extend at trace amounts farther out 
than this boundary, the directly observed a{r) curves (Figure 1) have plainly almost leveled 
off there, indicating that we are already including the main portion of the GCS. 
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Figure 2 shows log-log plots of the background-subtracted density profiles for the four 
BCGs, and Table 1 lists our deduced values for the power-law index b, as well as the adopted 
afjg. We note the bumpiness of the (Tci{r) curve for NGC 708 in Abell 262. It exhibits a clear 
peak around 1/25, and a less clear peak around 0/6. If not simply statistical fluctuations, 
these could be the result of uncompensated background contamination (dumpiness in the 
distribution of faint, ultra-distant galaxy clusters, or even some residual globular cluster 
populations which belong to neighboring galaxies in A262; there are several smaller elliptical 
galaxies as well as a companion spiral in the WFPC2 field) . 

Also, the slope of the (Tc;(r) curve for IC 4296 in A3565 declines smoothly out to about 
1.'5, and then exhibits a relatively sharp drop-off. This effect is not apparent in the other 
three systems. It is possible that this feature, if real (and not just a prosaic result of 
overestimating the background count), is due to tidal truncation. There are no large nearby 
galaxies in the immediate vicinity of IC 4296. However, Mulchaey et al. (1996) in a ROSAT 
study of groups and poor clusters find that the X-ray contours in A3565 arc not centered 
on the BCG, but rather at a point roughly halfway between IC 4296 and the neighboring 
spiral IC 4299. The mass that Mulchaey et al. (1996) have derived for the X-ray gas is 
about one-third of the mass of the entire cluster derived from the galaxy velocity distribution 
(Willmer et al. 1999). It might therefore be possible that the BCG oscillates about the 
center of the potential well of the cluster, which could be responsible for the truncation of 
its halo. 

For luminous giant E galaxies like these, the GCS profile shape is expected to be rather 
flat, becoming more shallow as galaxy luminosity increases. A general relation for the power- 
law slope is (Kaisler et al. 1996). 

b = -0.29 - 8.00 

The only one of our four target BCGs which deviates fairly strongly from this empirical 
relation is NGC 708 (an observed slope b — —1.0 versus an expected one of —1.7). It is, 
however, the only one in our sample with a cD envelope, which has clearly made it more 
radially extended than normal. 

3.2. Globuleir Cluster Luminosity Function 

The globular cluster luminosity function (GCLF) is conventionally defined as the number 
of globular clusters per unit magnitude. In this form it has a unimodal and nearly symmetric 
shape, and for analytical purposes can be reasonably approximated by a Gaussian function 
(cf. Harris (2001) for extensive discussion and background). For giant elliptical galaxies the 
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Fig. 1. — Radial profiles of the detected starfike objects around the four target BCGs. The 
dotted lines indicate the adopted background density levels, defined as the mean of the 
outermost annuli. 
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Fig. 2. — Log-log plot of radial profiles of the detected starlike objects around the four target 
BCGs. The background density, a^g, has been subtracted. 
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Table 1: Radial Density Profile Data 



Galaxy 




b 




(Jbg 




N708 (A262) 


-1.00 


± 


0.11 


11.6 ± 


3.1 


N5193 (A3560) 


-1.38 


± 


0.07 


2.8 ± 


1.9 


14296 (A3565) 


-1.02 


± 


0.08 


15.0 ± 


3.9 


N7014 (A3742) 


-1.71 


± 


0.11 


16.0 ± 


3.6 



Gaussian dispersion ay of the distribution is found to be 1.36 ± 0.03 (see below) while the 
peak frequency, also known as the GCLF turnover point, is equal to My — —7.33 ± 0.04 on 
a distance scale where the calibrating ellipticals (M87 and others in the Virgo cluster) are at 
d — 16 Mpc (see Harris (2001): Kavelaars et al. (2000)). The demonstrated small galaxy- 
to-galaxy spread in My makes it an interesting standard candle for distance measurements, 
because it is detectable to distances beyond which peculiar motions of galaxies bias the 
Hubble flow (see §3.5). 

To plot the GCLFs for our four GCSs, we have used completeness-corrected totals in 
0.25 mag bins above the 50% completeness level magnitudes. A background GCLF was sta- 
tistically subtracted, defined from the objects (mostly background galaxies) in the outskirts 
of the fields (same region as the one used to define a^g above). The best-fit Gaussian function 
was obtained with a constrained minimization method. 

When solving for the two important GCLF parameters (the apparent magnitude of the 
turnover, in this case and the dispersion o"/), it is important to note that their solutions 
are correlated (see Seeker & Harris (1993); Hanes & Whittaker (1987); Kavelaars et al. 
(2000)). If, as is true in our case, the photometric limit I{lim) is close to or just past the 
true turnover magnitude, the fitted values for the turnover and dispersion are constrained 
from only one side (the bright end). In this situation, attempts to solve simultaneously for 
the turnover and the dispersion tend to produce overestimates of both quantities. One can 
obtain a systematically more accurate fit by adopting a fiducial value for the dispersion, and 
solving only for the turnover magnitude. Fortunately, it has been shown (see Whitmore 
(1996) and Harris (2001)) for more than a dozen gE galaxies with well measured GCLFs, 
that the dispersion is very consistent from one galaxy to another, with (cy) = 1.36 ± 0.10 
rms scatter. We adopt that dispersion value here and solve only for the turnover magnitudes 
7°. 

Our full completeness-corrected and background-subtracted luminosity functions, along 
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with the fitted Gaussian functions (assuming ui — 1.36) are shown in Figure 3. The individ- 
ual points were weighted as {n/e{n)Y where n and e(n) denote the number of objects in the 
bin and the internal uncertainty. Table 2 lists the resulting turnover /-magnitudes obtained 
for the GCLFs, with successive columns giving the result from (a) only the PCI data, (b) 
only the WF data, (c) all data combined, and finally (d) the values determined previously 
by Lauer et al. The values result from the fitted /° once we add our adopted mean color 
index (y — 7)o = 1.1 ± 0.1. In the Table, the quoted 7° values have had Galactic extinction 
subtracted. The quoted uncertainties in the Table represent only the internal uncertainty of 
the fitting procedure. 

Our results for the four BCGs yield turnover magnitudes quite similar to those obtained 
by Lauer et al. (1998) for A262, A3742, but fainter than their values for A3560 and A3565. 
Lauer et al. employed the same basic fitting curve (a Gaussian with dispersion 1.4 magni- 
tudes), but there are two notable differences between their analysis and ours. First, we have 
used all the available WFPC2 data (not just the PCI data), more than doubling the total 
cluster populations used in the fit. Second, there are important differences of detail in the 
fitting procedure. We define an empirical background LF from the outer parts of the images; 
then subtract this background from the total number counts; and then fit the (residual) 
LF with the model Gaussian function. By contrast, Lauer et al. fit a model to the raw 
data which includes both the assumed Gaussian GCLF shape and a model background LF 
which rises smoothly and exponentially with magnitude. Either procedure depends for its 
accuracy on the correctness of the assumed background LF particularly at faint levels. Thus, 
for example, if their background model overestimated the true background at or beyond the 
true GCLF turnover, it would yield an /° which was artificially too bright. In cases such as 
these where the photometric limit is very close to the actual turnover, the true (external) 
uncertainty in is likely to be closer to ±0.3 magnitude (cf. the references cited above). 
However, primarily because of the much larger GCLF sample size in our data and our locally 
determined backgrounds, we believe our turnover determinations to be improvements over 
the previous work.^ 

In Figure 4 we show the GCLF from our data broken separately into the WF and PCI 
data (and with separate fitted Gaussians for both subsets). Notably, the fits from the PCI 
data alone yield consistently brighter turnovers than do the deeper WF data, indicating 
the importance of obtaining photometric limits at or beyond the actual turnover to avoid 



^We have compared the GCLFs from only our PCI data directly with those of Lauer et al. (see their 
Figure 8). After correction for the different bin sizes, we find that their actual GCLF datapoints and ours 
agree quite closely, to well within the internal errors of cither set. Our PCI data do not reach as faint as 
theirs, however, indicating that our detection threshold was more conservative. 
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systematic errors. The WF data have hmiting magnitudes similar to those of Lauer et al., 

and the turnover magnitudes we deduced from the WF data alone are generally closer to 
those deduced by Lauer et al., as are those from the total (PC+WF) dataset. These latter 
fits are certainly to be preferred over the PCI data alone. 

We note that the datapoints in the faintest bin for A262, A3560 and A3742 in Figure 3 
and in the WF Figure 4 fall clearly below the fitted GCLF curve. This does not appear to be 
due to an incorrect completeness correction at that level. The only other obvious possibility 
is that our empirical background has been overestimated in the last bin, but it is unclear 
whether that is the case without deeper data to draw on. 

Table 2: GCLF Turnover Magnitudes 



Galax}' 


/"(PCI) 


/"(WF) 


/"(total) 


/"(Lauer et al.) 


N708 (A262) 


24.66 


25.90 


25.74 ± 0.17 


25.85 ± 0.25 


N5193 (A3560) 


24.62 


25.81 


25.81 ± 0.18 


25.12 ± 0.25 


14296 (A3565) 


24.66 


25.72 


25.54 ± 0.16 


24.72 ± 0.25^ 


N7014 (A3742) 


25.58 


25.71 


25.66 ± 0.19 


25.77 ± 0.25 



^ The turnover magnitude for IC 4296 listed by Lauer et al. (1998) 
in their Table 2 (/q = 25.72) is a typographical error. It should read 
7° = 24.72, which is the value from their graph (Figure 8) (Lauer 
2001). 



3.3. Luminosity and Mass Distribution Functions 

The luminosity distribution function (LDF), or number of clusters per unit luminosity 
{dN/dL), is the visible signature of the cluster mass distribution and is a more physically 
oriented representation of the GCLF. In this plane it usually appears as a rough power law, 
dN/ dL ~ L~°'^ . By adopting a mass-to-light ratio one can also immediately obtain the mass 
distribution function dN/dM ~ M~°"^ , where — oim as long as the M/L ratio is roughly 
independent of cluster mass (see Mandushev et al. (1991) and McLaughhn (2000)). 

It has been noted many times that within observational uncertainty, the globular cluster 
mass distribution follows a simple power law relation which has the same shape as the mass 
distribution of giant molecular clouds (GMCs) in large spirals, cloud cores embedded in 
GMCs, and giant HII regions in large spirals (e.g. Harris & Pudritz (1994) and subsequent 
papers). The power-law slope um falls consistently in the range 1.8 ± 0.2 for L > IO^Lq, a 
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Fig. 3. — The data defining the globular cluster luminosity function (GCLF) for bins brighter 
than the 50% completeness level. The solid fines are Gaussians of width cry = 1.36. 
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Fig. 4. — The data defining the globular cluster luminosity functions (GCLFs) for bins above 
the 50% completeness level on the PCI (dashed) and WF (solid) chips. The fitted lines are 
Gaussians of width cry = 1.36. 
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value which can be at least approximately explained by cluster formation pictures ranging 
from collisional growth of protocluster gas clouds (McLaughlin & Pudritz 1996) to turbulence 
spectra (Elmegreen & Efremov 1997). 

Here we define the LDFs observationally in a separate procedure from the GCLF: we 
convert the apparent magnitudes of the detected objects to absolute magnitudes My us- 
ing the galaxy distance moduli and adopted cluster color, and then to luminosities Ly. The 
background LDF is subtracted and the data binned in steps of IO^Lq. For the mass distribu- 
tion function we have multiplied the luminosity values by a mass-to-light ratio M/ Ly = 1.45 
(McLaughlin 2000) and binned the data in steps of IO^Mq. 

We note that the lowest luminosity (and hence mass) point for each of the functions has 
been excluded from the fits. This corresponds to a point L < lO^L©, since at this luminosity, 
the slope of the power law relation changes for the distribution function (see McLaughlin 
(1994)). 

Figure 5 shows the LDFs which we have obtained for the four galaxies. Table 3 lists 
the values obtained for the slopes of the power law relations of the luminosity and mass 
distribution functions, ai and aM- (These are slightly but not significantly different in each 
pair because of the different bin boundary locations in L and M.) With the exception of 
NGC 5193 in A3560 which has a value of ckm notably lower than the expected 1.8 ± 0.2, the 
other three galaxies are entirely consistent with the other giant ellipticals in Virgo, Fornax 
and elsewhere. 

Table 3: Luminosity and Mass Distribution Function Data 



Galaxy 


OIL 


aM 


N708 (A262) 


-1.97 ± 0.11 


-1.86 ± 0.08 


N5193 (A3560) 


-1.57 ± 0.07 


-1.48 ± 0.05 


14296 (A3565) 


-1.85 ± 0.07 


-1.74 ± 0.07 


N7011 (A3712) 


-1.96 ± 0.10 


-1.91 ± 0.08 



3.4. Total Populations and Specific Frequencies 

Next we estimate the total populations and specific frequencies of the GCSs. Here Sj^ 
is the number of clusters per unit galaxy luminosity (Harris & van den Bergh 1981; Harris 
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Fig. 5. — The data defining the luminosity distribution function (LDF). The sohd hnes are 
the least squares fits, with slopes given in Table 3. The lowest luminosity point has been 
excluded from the fit in all cases (see text). 
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et al. 1991): 

5^ = 7V^,xlO°-^(^?+i5) 

where is the total number of clusters, and My is the integrated absolute magnitude of 
the host galaxy. Sn is found to differ among gE galaxies by more than an order of magnitude 
(Harris, Harris & McLaughlin 1998; Harris 2001). It is known that for cD-typc and BCG 
galaxies the specific frequency increases systematically with the total galaxy luminosity, the 
size of the surrounding cluster of galaxies, and the X-ray halo gas mass (Blakeslee et al. 
1997; Blakeslee 1999; Harris, Harris & McLaughlin 1998; McLaughlin 1999; Kavelaars 
1999), albeit with a factor-of-two scatter that is not understood in detail. 

We calculate the total number of clusters by integrating the radial profiles and scaling 
up the result by the fraction of the area under the bright half of the GCLFs (the complete- 
ness magnitude used for the radial profiles corresponds to 90% completeness, while for the 
GCLFs it is 50%). Following previously established convention (Harris 2001), the result 
is then doubled, which imphcitly assumes that the GCLF is symmetric about the turnover 
magnitude. We can therefore think of specific frequency as equivalent to the number of 
bright clusters in a galaxy. Also, the value of is fairly insensitive to the assumptions in 
the galaxy distance, because changes in distance will affect the calculated galaxy himinosity 
and total cluster population in the same sense (see Harris & van den Bergh (1981)). 

The total populations along with the specific frequencies are listed in Table 4. Fig- 
ure 6 shows our results along with those from other gE galaxies previously published. 

NGC 708 in A262 has the highest value for Sn, consistent with the fact that it is the 
only genuine cD-type galaxy of the four studied. The other three galaxies are central BCGs 
without the extended cD envelope, and have lower values for Sn- Somewhat of an outlier is IC 
4296 in A3565, with a Sn value of 2.6, definitely on the low end for a such a high-luminosity 
elliptical. Its Sn is, however, similar to those of "field" ellipticals which are commonly 
thought to have formed by major mergers between relatively less cluster-rich disk galaxies. 
Even if the merging galaxies are quite gas-rich, a high— S'at elliptical would not necessarily 
result, since new field stars and star clusters both form during the merger, and the net ratio 
of clusters to field stars in the final merger product could either increase or decrease. The 
final Sn would be higher only if the efficiency of cluster formation was considerably enhanced 
over the cluster formation that took place in the original protogalactic epoch. In actually 
observed cases of recent disk/disk mergers, what appears to be emerging in every case is an 
elliptical with S'at ~ 2 (see Harris (2001) for more extensive discussion). 
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Fig. 6. — Specific Frequency Sjsf plotted against luminosity for elliptical galaxies. Squares 
represent our results for the BCGs, while open circles are results for other BCGs and gEs. 
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Table 4: Specific Frequency Data 



Galaxy 


Mv 


Nci 




N708 (A262) 


-22.14 


5924 ± 302 


8.24 ± 1.65 


N5193 (A3560) 


-22.75 


5878 ± 355 


4.66 ± 0.93 


14296 (A3565) 


-23.49 


6400 ± 265 


2.58 ± 0.52 


N7014 (A3742) 


-21.99 


1634 ± 88 


2.62 ± 0.52 



3.5. The Hubble Diagram 

By using our results for the GCLF turnover magnitudes, along with previous results 
for several other gE galaxies from the literature a classic "Hubble Diagram" of redshift cz 
against the apparent turnover magnitude of the GCLF can be constructed (see Harris 
(2001) and Kavelaars et al. (2000) for the first use of this diagram for GCLFs). 

Hubble's law states cz — Hod. This can be rewritten as: 

log(c^) = 0.2V° + log//o - 0.2M^ - 5 

where My is the GCLF turnover luminosity for gE galaxies, and Hq is expressed in the 
usual units of km s~^ Mpc~^. Plotting \og{cz) against the apparent magnitude gives a 
straight line of slope 0.2, and a zcropoint which contains Hq and My. 

The available data for a total of eleven BCG galaxies or groups are listed in Table 5 
and plotted in Figure 7. The values for the Virgo and Fornax clusters arc the weighted 
mean (V^) values of the individual galaxies (sec Kavelaars ct al. (2000)). The mean radial 
velocities of each galaxy or group are taken from Fabcr ct al. (1989), Girardi et al. (1993), 
Huchra (1988), Binggeli, Popescu & Tammann (1993), Hamuy et al. (1996), CoUess & 
Dunn (1996), and Lauer et al. (1998). Also, the recession velocities cz for the target galaxies 
assume a Local Group infall to Virgo of 250 ± 100 km s"^ (e.g.. Ford et al. (1996); Hamuy 
et al. (1996); Jerjen & Tammann (1993), among many others). 

The four galaxies we have investigated here fall well within the pattern established by 
the others (from Virgo at low redshift out to Coma at the highest redshift) . A value of Hq 
near ~ 70 km Mpc~^, which represents a consensus of recent determinations (Preedman 
et al. 2001), along with = -7.33 (Kavelaars et al. 2000; Harris 2001), matches the 
total range of points within their measurement uncertainties. 
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22 24 26 28 



GCLF Turnover Magnitude V° 

Fig. 7. — Hubble diagram for GCLF turnover magnitudes. The galaxy redshift cz is plotted 
against the apparent magnitude of the GCLF turnover, . The solid line is the expected 
relation for Hq = 70 and an absolute turnover magnitude My = —7.33 (see text). 
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Table 5: GCLF turnover magnitudes in distant ellipticals 



Cluster 


Galaxy 


Redshift cz (km/s) 


(GCLF) 


Source 


Virgo 


6 gEs 


1300 


23.71 ± 0.03 


1 


Fornax 


6 gEs 


1400 


23.85 ±0.04 


1 


NGC 5846 


NGC 5846 


2300 


25.08 ±0.10 


2 


Coma 


IC 4051 


7100 


27.75 ±0.20 


3 


Coma 


NGC 4874 


7100 


27.82 ±0.12 


1 


Coma 


NGC 4881 


7100 


> 27.6 


4 


Coma 


NGC 4926 


7100 


27.90 ±0.20 


5 


A 262 


NGC 708 


4650 


26.84 ±0.17 


6 


A 3560 


NGC 5193 


4020 


26.91 ±0.18 


6 


A 3565 


IC 4296 


4110 


26.63 ±0.16 


6 


A 3742 


NGC 7014 


4680 


26.75 ±0.19 


6 



Sources: (1) Kavelaars et al. (2000); (2) Forbes (1996); (3) Woodworth & Harris (2000); 
(4) Baum et al. (1995); (5) Kavelaars (2001); (6) this paper 

4. Summciry 

Using deep /-band photometry from the /JS'T/WFPC2 archive we have studied the 
globular cluster systems in NGC 708, NGC 5193, IC 4296 and NGC 7014, the BCGs in 
A262, A3560, A3565 and A3742. The photometry allowed us to construct the globular 
cluster luminosity functions for each of the above cluster systems, reaching to or just past 
the turnover point. 

We have used a constrained method to fit a Gaussian function to each of the four 
GCLFs, adopting a width uy = 1.36. The resulting background and extinction corrected 
y-band turnover magnitudes were found to be at = 26.84 ± 0.17, 26.91 ± 0.18, 26.63 
± 0.16, 26.75 ± 0.19 respectively. These values improve on previous results of Lauer et al. 
(1998). The Hubble diagram generated from our GCLF turnover data combined with other 
material from the literature matches Hq ~ 70 km s~^ Mpc~^. 

We have obtained the luminosity and mass distribution functions for the GCSs. The 
slopes of the power law relations were found to be consistent with those for other giant 
elliptical galaxies, with the exception of NGC 5193 in A3560. 

The total cluster populations and the specific frequencies Sn for each system were 
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calculated. The specific frequency of IC 4296 in A2565 is surprisingly low for its high 
luminosity, while the values for the other three galaxies fall within the established trend of 
specific frequency versus galaxy luminosity 

We would hke to thank J.J.Kavelaars and Marcel VanDalfsen for useful comments and 
advice. This work was supported by the Natural Sciences and Engineering Research Council 
of Canada through a research grant to WEH. 
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